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1. Introduction 

The 3+1 formalism is the basis of most modern numerical relativity and has lead, 
along with alternative approaches [82] , to the recent successes in the binary black hole 
merger problem [7J ESI E3 E3 EH EE] (see El] [69] EB] for a review). Thanks to 
the 3+1 formalism, the resolution of Einstein equation amounts to solving a Cauchy 
problem, namely to evolve "forward in time" some initial data. However this is a 
Cauchy problem with constraints. This makes the set up of initial data a non trivial 
task, because these data must fulfill the constraints. In this lecture, we present the 
most wide spread methods to deal with this problem. Notice that we do not discuss 
the numerical techniques employed to solve the constraints (see e.g. Choptuik's lecture 
for finite differences [32] and Grandclement and Novak's review for spectral methods 
[58]). 

Standard reviews about the initial data problem are the articles by York [106] and 
Choquet-Bruhat and York [36 a . Recent reviews are the articles by Cook [37], Pfeiffer 
[TP] and Bartnik and Isenberg [ID] . 

2. The initial data problem 

2.1. 3+1 decomposition of Einstein equation 

In this lecture, we consider a spacetime (M,g), where M. is a four-dimensional 
smooth manifold and g a Lorentzian metric on M. . We assume that (M , g) is globally 
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hyperbolic, i.e. that Ai can be foliated by a family (E 4 ) t6 R of spacelike hypersurfaces. 
We denote by 7 the (Riemannian) metric induced by g on each hypersurface E t and 
K the extrinsic curvature of E t , with the same sign convention as that used in the 
numerical relativity community, i.e. for any pair of vector fields (u, v) tangent to E t , 
g(u, V^n) = —K(u,v), where n is the future directed unit normal to E t and V is 
the Levi-Civita connection associated with g. 

The 3+1 decomposition of Einstein equation with respect to the foliation (St)tgR 
leads to three sets of equations: (i) the evolution equations of the Cauchy problem 
(full projection of Einstein equation onto St), (ii) the Hamiltonian constraint (full 
projection of Einstein equation along the normal n), (iii) the momentum constraint 
(mixed projection: once onto E t , once along n). The latter two sets of equations do 
not contain any second derivative of the metric with respect to t. They are written^] 

R + K 2 - KijK ij = IGttE (Hamiltonian constraint), (1) 
D J Kij — DiK — %irpi (momentum constraint), (2) 

where R is the Ricci scalar (also called scalar curvature) associated with the 3-metric 
7, K is the trace of K with respect to 7: K = 7 y \Kjj, D stands for the Levi-Civita 
connection associated with the 3-metric 7, and E and pi are respectively the energy 
density and linear momentum of matter, both measured by the observer of 4-velocity 
n (Eulerian observer). In terms of the matter energy-momentum tensor T they are 
expressed as 

E = T ilv n»n v and Pi = -T^n^i- (3) 

Notice that Eqs. U)-© involve a single hypersurface Eg, not a foliation (E t ) t R . In 
particular, neither the lapse function nor the shift vector appear in these equations. 

2.2. Constructing initial data 

In order to get valid initial data for the Cauchy problem, one must find solutions to 
the constraints {1} and ©. Actually one may distinguish two problems: 

• The mathematical problem: given some hypersurface So, find a Riemannian 
metric 7, a symmetric bilinear form K and some matter distribution (E, p) on So 
such that the Hamiltonian constraint |T]) and the momentum constraint © are 
satisfied. In addition, the matter distribution (E, p) may have some constraints 
from its own. We shall not discuss them here. 

• The astrophysical problem: make sure that the solution to the constraint equations 
has something to do with the physical system that one wish to study. 

Facing the constraint equations ([1]) and ([2]), a naive way to proceed would be to 
choose freely the metric 7, thereby fixing the connection D and the scalar curvature 
R, and to solve Eqs. CO)-© for K- Indeed, for fixed 7, E, and p, Eqs. U)-© form 
a quasi-linear system of first order for the components Kij. However, as discussed 
by Choquet-Bruhat [45] , this approach is not satisfactory because we have only four 
equations for six unknowns Kij and there is no natural prescription for choosing 
arbitrarily two among the six components Kij . 

In 1944, Lichnerowicz [70] has shown that a much more satisfactory split of 
the initial data (7, K) between freely choosable parts and parts obtained by solving 

I we are using the standard convention for indices, namely Greek indices run in {0, 1, 2, 3}, whereas 
Latin ones run in {1, 2, 3} 
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Eqs. (HJ-© is provided by a conformal decomposition of the metric 7. Lichnerowicz 
method has been extended by Choquet-Bruhat (1956, 1971) [351 133], by York and 
6 Murchadha (1972, 1974, 1979) PM QUI US [TUB] and more recently by York and 
Pfeiffer (1999, 2003) |107U80j . Actually, conformal decompositions are by far the most 
widely spread techniques to get initial data for the 3+1 Cauchy problem. Alternative 
methods exist, such as the quasi-spherical ansatz introduced by Bartnik in 1993 [8] or 
a procedure developed by Corvino (2000) [39] and by Isenberg, Mazzeo and Pollack 
(2002) [63] for gluing together known solutions of the constraints, thereby producing 
new ones. Here we shall limit ourselves to the conformal methods. 



2. 3. Conformal decomposition of the constraints 

In the conformal approach initiated by Lichnerowicz [70] . one introduces a conformal 
metric 7 and a conformal factor ^> such that the (physical) metric 7 induced by the 
spacetime metric on the hypersurface St is 

7<i=* 4 7tf- (4) 

We could fix some degree of freedom by demanding that detTy = 1. This would 
imply \& = (det 7i J ) 1 / 12 . However, in this case 7 and \& would be tensor densities. 
Moreover the condition det 7y = 1 has a meaning only for Cartesian-like coordinates. 
In order to deal with tensor fields and to allow for any type of coordinates, we proceed 
differently and introduce a background Riemannian metric / on S t . If the topology of 
St allows it, we shall demand that / is flat. Then we replace the condition det^ = 1 
by det 7y = det /y . This fixes 

*=fdet^y /12 - 



Kij =: Aij + -Kjij. (6) 



^det/y 

\I> is then a genuine scalar field on S t (as a quotient of two determinants) . Consequently 
7 is a tensor field and not a tensor density. 

Associated with the above conformal transformation, there are two decomposi- 
tions of the traceless part of the extrinsic curvature, the latter being defined by 

1 

3~ 

These two decompositions are 

.1'-' : >l< : ".K (7) 

A ij : >l< -A'*. (8) 

The choice —10 for the exponent of \& in Eq. ([7]) is motivated by the following identity, 
valid for any symmetric and traceless tensor field, 

DjA i3 = y- 10 Dj (V w A ij ) , (9) 

where Dj denotes the covariant derivative associated with the conformal metric 7. 
This choice is well adapted to the momentum constraint, because the latter involves 
the divergence of K. The alternative choice, i.e. Eq. ©, is motivated by time 
evolution considerations, as we shall discuss below. For the time being, we limit 
ourselves to the decomposition ([7]), having in mind to simplify the writing of the 
momentum constraint. 
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By means of the decompositions (QJ, ^ and ([7]), the Hamiltonian constraint (fT]) 
and the momentum constraint © are rewritten as (see Ref. [FT] for details) 

Di&V - -m + -A i0 A i3 ty- 7 + 2nE^- 3 - —K 2 V° = 0, (10) 
8 8 12 

DjA* '- -^ 6 &K = 87TP 4 , (11) 
3 

where i? is the Ricci scalar associated with the conformal metric 7 and we have 
introduced the rescaled matter quantities 

E := and p l := ^ w p\ (12) 

Equation (|10[) is known as Lichnerowicz equation, or sometimes Lichnerowicz- York 
equation. The definition of p l is such that there is no \& factor in the right-hand side 
of Eq. (TTT]) . On the contrary the power 8 in the definition of E is not the only possible 
choice. As we shall see in § 13.41 it is chosen (i) to guarantee a negative power of ^> in 
the E term in Eq. (fT0|) , resulting in some uniqueness property of the solution and (ii) 
to allow for an easy implementation of the dominant energy condition. 

3. Conformal transverse-traceless method 

3.1. Longitudinal/transverse decomposition of A 13 

In order to solve the system (TlU])-([ni), York (1973,1979) [TTil[Tn51[M] has decomposed 
A 13 into a longitudinal part and a transverse one, setting 

A ij = (Lxyi +A% T , (13) 

where Al^ T is both traceless and transverse (i.e. divergence-free) with respect to the 
metric 7: 

j i:j A^ T = and DjAfi T = 0, (14) 

and (LXy j is the conformal Killing operator associated with the metric 7 and acting 
on the vector field X: 

(Lxy 3 := D i X j + b 3 x r - -D k X k f 3 . (15) 

3 

(LXyi is by construction traceless: 

%{LX) i3 =0 (16) 

(it must be so because in Eq. (fT3|) both A lJ and Al£, T are traceless). The kernel of 
L is made of the conformal Killing vectors of the metric 7, i.e. the generators of 
the conformal isometries (see e.g. Ref. [51] for more details). The symmetric tensor 
(LX) lJ is called the longitudinal part of A* 3 , whereas Al^ T is called the transverse part. 

Given A 13 , the vector X is determined by taking the divergence of Eq. (fl3)) : 
taking into account property (fl~4|) . we get 

Dj(LXy j =DjA i3 . (17) 

The second order operator Dj(LXy 3 acting on the vector X is the conformal vector 
Laplacian A^: 

A L X 1 := Dj(LX) i: > = D J D j X i + -DWjX 3 + R^X 3 , (18) 
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where the second equality follows from the Ricci identity applied to the connection _D, 
Rij being the associated Ricci tensor. The operator Ax, is elliptic and its kernel is, in 
practice, reduced to the conformal Killing vectors of 7, if any. We rewrite Eq. (fT7)l as 
A L X l = DjAv. (19) 

The existence and uniqueness of the longitudinal/transverse decomposition (|13[) 
depend on the existence and uniqueness of solutions X to Eq. (fT9|) . We shall consider 
two cases: 

• Eo is a closed manifold, i.e. is compact without boundary; 

• (£0,7) is an asymptotically flat manifold, i.e. is such that the background metric 
/ is flat (except possibly on a compact sub-domain B of E t ) and there exists a 
coordinate system (x l ) — (x,y,z) on £ t such that outside B, the components 
of / are /y = diag(l,l,l) ("Cartesian-type coordinates") and the variable 
r := y/ x 2 + y 2 + z 2 can take arbitrarily large values on E t ; then when r — ► +00, 
the components of 7 and K with respect to the coordinates (x l ) satisfy 

7tf = /«+0(r- 1 ) and ^ = 0(r" 2 ), (20) 

Kii = 0(r- 2 ) and ^ = ( r -3). (2 1) 

In the case of a closed manifold, one can show (see Appendix B of Ref. [SI] for details) 
that solutions to Eq. (fT9"|) exist provided that the source DjA 1 ^ is orthogonal to all 
conformal Killing vectors of 7, in the sense that 

VCekcrZ, f j ij C i D k A :ik ^d 3 x = Q. (22) 

But the above property is easy to verify: using the fact that the source is a pure 
divergence and that Eo is closed, we may integrate the left-hand side by parts and 
get, for any vector field C, 



%C l D k A^y/^fd d x = -- j ij % l (LCy k A^^d i x. (23) 

Then, obviously, when C is a conformal Killing vector, the right-hand side of the 
above equation vanishes. So there exists a solution to Eq. (fT9)) and this solution is 
unique up to the addition of a conformal Killing vector. However, given a solution 
X, for any conformal Killing vector C, the solution X + C yields to the same value 
of LX, since C is by definition in the kernel of L. Therefore we conclude that the 
decomposition (fl3|) of A % i is unique, although the vector X may not be if (£0,7) 
admits some conformal isometries. 

In the case of an asymptotically flat manifold, the existence and uniqueness is 
guaranteed by a theorem proved by Cantor in 1979 [3D] (see also Appendix B of 
Ref. [87] as well as Refs. [35] [51]). This theorem requires the decay condition 

Sh=°^ < 24 > 

in addition to the asymptotic flatness conditions (|2TH) . This guarantees that 

^■=0(r- 3 ). (25) 
Then all conditions are fulfilled to conclude that Eq. (TT9|) admits a unique solution X 
which vanishes at infinity. 

To summarize, for all considered cases (asymptotic flatness and closed manifold) , 
any symmetric and traceless tensor A^ (decaying as 0(r~ 2 ) in the asymptotically flat 
case) admits a unique longitudinal/transverse decomposition of the form (fT3| . 
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3.2. Conformal transverse-traceless form of the constraints 

Inserting the longitudinal/transverse decomposition Q13[) into the constraint equations 
pop and (fTTj) and making use of Eq. (fH)|) yields to the system 

1 ^ 1 



d.d 1 ^ - -m + - 

8 8 



1 



-7 



where 



+ 2ttE^- 3 - — i^ 2 * 5 = 0, (26) 

A L X 1 - -^ 6 D l K = 8np\ (27) 
3 

(LX^-.^^^LXf 1 and Aj? := W^tt- (28) 

With the constraint equations written as (f2"6")l and (|27p. we see clearly which part 
of the initial data on So can be freely chosen and which part is "constrained" : 

• free data: 

— conformal metric 7; 

— symmetric traceless and transverse tensor A^ T (traceless and transverse are 
meant with respect to 7: "fijAl^ T — and DjA^ T = 0); 

— scalar field K; 

— conformal matter variables: (E,p l ); 

• constrained data (or "determined data"): 

— conformal factor obeying the non-linear elliptic equation (|26|) 
(Lichnerowicz equation) 

— vector X, obeying the linear elliptic equation (j27|) . 

Accordingly the general strategy to get valid initial data for the Cauchy problem is 
to choose (jij, A l r£ T , K, E,p l ) on Eo and solve the system (|26|) - ([27| to get VP and X 1 . 
Then one constructs 



Hi 






3 


(29) 


K ij 


= *- 10 {[ixy 1 - 




(30) 


E 


= ^-*E 






(31) 


P* 


= *- 10 p* 






(32) 



and obtains a set (-y,K,E,p) which satisfies the constraint equations ([I])-©. This 
method has been proposed by York (1979) [106] and is naturally called the conformal 
transverse traceless (CTT) method. 

3. 3. Decoupling on hypersurfaces of constant mean curvature 

Equations (|26p and (f2"7) are coupled, but we notice that if, among the free data, we 
choose if to be a constant field on Eo, 

K = const, (33) 

then they decouple partially : condition (f3"3")) implies D l K = 0, so that the momentum 
constraint (|2T|) becomes independent of 

A L X l = 8ttp 1 (K = const). (34) 
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The condition (|33|) on the extrinsic curvature of So defines what is called a constant 
mean curvature (CMC) hypersurface. Indeed let us recall that K is nothing but 
(minus three times) the mean curvature of (£ ,~/) embedded in (A4,g). A maximal 
hypersurface, having K — 0, is of course a special case of a CMC hypersurface. On a 
CMC hypersurface, the task of obtaining initial data is greatly simplified: one has first 
to solve the linear elliptic equation (|34ll to get X and plug the solution into Eq. ([26]) 
to form an equation for 4". Equation (|34p is the conformal vector Poisson equation 
discussed above (Eq. (fT9|) . with DjA 1 ? replaced by 8717?). We know then that it is 
solvable for the two cases of interest mentioned in Sec. 13.11 closed or asymptotically 
flat manifold. Moreover, the solutions X are such that the value of LX is unique. 

3.4- Lichnerowicz equation 

Taking into account the CMC decoupling, the difficult problem is to solve Eq. ([26]) 
for ty. This equation is elliptic and highly non-linea4§|. It has been first studied 
by Lichnerowicz [70] [7T] in the case K = (So maximal) and E = (vacuum). 
Lichnerowicz has shown that given the value of \1/ at the boundary of a bounded domain 
of S (Dirichlet problem), there exists at most one solution to Eq. ([2l>]) . Besides, he 
showed the existence of a solution provided that A^A % ^ is not too large. These early 
results have been much improved since then. In particular Cantor [29j has shown that 
in the asymptotically flat case, still with K = and E = 0, Eq. (|26p is solvable if 
and only if the metric 7 is conformal to a metric with vanishing scalar curvature (one 
says then that 7 belongs to the positive Yamabe class) (see also Ref. [74]). In the 
case of closed manifolds, the complete analysis of the CMC case has been achieved by 
Isenberg (1995) [62] , 

For more details and further references, we recommend the review articles by 
Choquet-Bruhat and York 36J and Bartnik and Isenberg [TO]. Here we shall simply 
repeat the argument of York |107j to justify the rescaling (fl"2]) of E. This rescaling is 
indeed related to the uniqueness of solutions to the Lichnerowicz equation. Consider 
a solution \& to Eq. ([2T>]) in the case K = 0, to which we restrict ourselves. Another 
solution close to \l/o can be written \& = "Jo + e , with |e| <C 

Di&^o + e) - ^(* + e) + \h,A^ (* + e)" 7 + 27r£(* + e)~ 3 = 0. (35) 

o o 

Expanding to the first order in e/^o leads to the following linear equation for e: 

Dib l e~ ae = 0, (36) 

with 

a := ±R + gi«l«>o 8 + 6^*0 4 - ( 37 ) 

Now, if a > 0, one can show, by means of the maximum principle, that the solution 
of ([36]) which vanishes at spatial infinity is necessarily e = (see Ref. [34] or § B.l of 
Ref. [3S]). We therefore conclude that the solution ^ to Eq. (|2"6"]) is unique (at least 
locally) in this case. On the contrary, if a < 0, non trivial oscillatory solutions of 
Eq. ([36]) exist, making the solution ^0 not unique. The key point is that the scaling 
(|12|1 of E yields the term +67r_E5' ~ 4 in Eq. ([37]) . which contributes to make a positive. 
If we had not rescaled E, i.e. had considered the original Hamiltonian constraint, the 

§ although it is quasi-linear in the technical sense, i.e. linear with respect to the highest-order 
derivatives 
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contribution to a would have been instead — IQitE^q, i.e. would have been negative. 
Actually, any rescaling E = ^5> S E with s > 5 would have work to make a positive. The 
choice s = 8 in Eq. (fTZ|) is motivated by the fact that if the conformal data (E^p 1 ) 
obey the "conformal" dominant energy condition 



E>^vP l p3, (38) 

then, via the scaling (TT2]) of <p % , the reconstructed physical data (E,p l ) will 
automatically obey the dominant energy condition 



/'• ' \ '.jl'P'- (39) 

4. Conformally flat initial data by the CTT method 

4-1- Momentarily static initial data 

In this section we search for asymptotically flat initial data (£o,7,-f0 by the CTT 
method exposed above. As a purpose of illustration, we shall start by the simplest 
case one may think of, namely choose the freely specifiable data (7^, A^ Tl K, E^p 1 ) 
to be a flat metric: 

lij=fij> (40) 
a vanishing transverse-traceless part of the extrinsic curvature: 

A% T = 0, (41) 
a vanishing mean curvature (maximal hypersurface) 

K = 0, (42) 
and a vacuum spacetime: 

£7 = 0, p l = 0. (43) 

Then Di = T>i, where T> denotes the Levi-Civita connection associated with f,R = 
(f is flat) and the constraint equations ([26]) - ([27]) reduce to 

A* + I {LX)ij {EX) 13 = (44) 
8 

A L X l = 0, (45) 

where A and A^ are respectively the scalar Laplacian and the conformal vector 
Laplacian associated with the flat metric /: 

A:=V{D l and A L X l := VjV 3 X 1 + ^-V'VjX 1 . (46) 

Equations ([4~H) - P5)l must be solved with the boundary conditions 

= 1 when r — > 00 (47) 

X = when r — > 00, (48) 

which follow from the asymptotic flatness requirement. The solution depends on the 
topology of So, since the latter may introduce some inner boundary conditions in 
addition to (|4"7|) - (|4"5)l . 

Let us start with the simplest case: So = K 3 . Then the unique solution of Eq. (|45)l 
subject to the boundary condition (|48|) is 

X = 0. (49) 
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Figure 1. Hypcrsurface Eo as R 3 minus a ball, displayed via an embedding 
diagram based on the metric -y, which coincides with the Euclidean metric on 
R 3 . Hence So appears to be fiat. The unit normal of the inner boundary S with 
respect to the metric ~f is s. Notice that D ■ s > 0. 



Consequently (LX)^ — 0, so that Eq. (14"4")1 reduces to Laplace equation for 4": 

A* = 0. (50) 

With the boundary condition ([171) . there is a unique regular solution on R 3 : 

* = 1. (51) 

The initial data reconstructed from Eqs. ([29]) -([30 |) is then 

7 = / (52) 
K = 0. (53) 

These data correspond to a spacelike hyperplane of Minkowski spacetime. 
Geometrically the condition K = is that of a totally geodesic hyper surf ace [i.e. all 
the geodesies of (E t ,7) are geodesies of (M,g)]. Physically data with K = are said 
to be momentarily static or time symmetric. Indeed, if we consider a foliation with 
unit lapse around Eo (geodesic slicing), the following relation holds: C n g = —2K, 
where C n denotes the Lie derivative along the unit normal n. So if K = 0, C n g = 0. 
This means that, locally (i.e. on Eo), n is a spacetime Killing vector. This vector 
being timelike, the configuration is then stationary. Moreover, the Killing vector n 
being orthogonal to some hypersurface (i.e. Eo), the stationary configuration is called 
static. Of course, this staticity properties holds a priori only on Eo since there is no 
guarantee that the time development of Cauchy data with K — at t = maintains 
K = at t > 0. Hence the qualifier 'momentarily' in the expression 'momentarily 
static' for data with K = 0. 



4-2. Slice of Schwarzschild spacetime 

To get something less trivial than a slice of Minkowski spacetime, let us consider a 
slightly more complicated topology for Eo, namely M 3 minus a ball (cf. Fig. [T]). The 
sphere S delimiting the ball is then the inner boundary of Eo and we must provide 
boundary conditions for \I> and X on S to solve Eqs. (j44])-(|45l). For simplicity, let us 
choose 

X\ s = 0. (54) 

Altogether with the outer boundary condition (|48|) . this leads to X being identically 
zero as the unique solution of Eq. (|45|) . So, again, the Hamiltonian constraint reduces 
to Laplace equation 

A* = 0. (55) 
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Figure 2. Same hypersurface So as in Fig. [T] but displayed via an embedding 
diagram based on the metric *y instead of •y. The unit normal of the inner 
boundary S with respect to that metric is s. Notice that D ■ s = 0, which 
means that S is a minimal surface of (£rj, T)- 



If we choose the boundary condition = 1, then the unique solution is ^ = 1 
and we are back to the previous example (slice of Minkowski spacetime). In order to 
have something non trivial, i.e. to ensure that the metric 7 will not be flat, let us 
demand that 7 admits a closed minimal surface, that we will choose to be S. This 
will necessarily translate as a boundary condition for W since all the information on 
the metric is encoded in ^ (let us recall that from the choice (j40|) . 7 = 4' 4 /). S is a 
minimal surface of (£0,7) iff its mean curvature vanishes, or equivalently if its unit 
normal s is divergence- free (cf. Fig. [2]): 

As*| s = 0. (56) 

This is the analog of V • n = for maximal hypersurfaces, the change from minimal 
to maximal being due to the change of metric signature, from the Riemannian to the 
Lorentzian one. Expressed in term of the connection D = T> (recall that in the present 
case 7 = /), condition (|56| is equivalent to 

V i {^ 6 s% = 0. (57) 

Let us rewrite this expression in terms of the unit vector s normal to S with respect 
to the metric 7 (cf. Fig. [I]); we have 

s = ^- 2 s, (58) 

since 7(s, s) = 5 ,_4 7(s, s) = ~f(s, s) = 1. Thus Eq. (|57p becomes 



4 s< 



= 0. (59) 

s 

Let us introduce on So a coordinate system of spherical type, (x l ) = (r,9,ip), such 
that (i) fij = diag(l, r 2 , r 2 sin 2 9) and (ii) iS is the sphere r = a, where a is some 
positive constant. Since in these coordinates y/J = r 2 sin# and s l — (1,0,0), the 
minimal surface condition (|59p is written as 

il:(^ 2 ) =o, (60) 
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Figure 3. Extended hypersurface Eg obtained by gluing a copy of So at the 
minimal surface S\ it defines an Einstein- Rosen bridge between two asymptotically 
fiat regions. 



i.e. 

a7 + 2^ 



= (61) 



This is a boundary condition of mixed Newmann/Dirichlet type for '5 . The unique 
solution of the Laplace equation (f55|) which satisfies boundary conditions (|47|) and 

m is 

* = 1 + -. (62) 

r 

The parameter a is then easily related to the ADM mass m of the hypersurface So- 
Indeed for a conformally flat 3-metric (and more generally in the quasi-isotropic gauge, 
cf. Chap. 7 of Ref. [FT]), the ADM mass to is given by the flux of the gradient of the 
conformal factor at spatial infinity: 

1 / 

m= lim (b ——r • sin 6 d6 dip 

2tt r^oo J r=const dr 

= - — lim 47rr 2 — 1 + - = 2a. 63 
2tt i — >oo or V r / 

Hence a = raj 2 and we may write 

* = (64) 

2r 

Therefore, in terms of the coordinates (r, 6*, </?), the obtained initial data (7, K) are 



to \ 



4 



Jij = [I + y t ) diag(l,^,^sin0) (65) 

K tj = 0. (66) 

So, as above, the initial data are momentarily static. Actually, we recognize on (|65|) - 
(|66[) a slice i = const of Schwarzschild spacetime in isotropic coordinates. 

The isotropic coordinates (r, 0, </?) covering the manifold S are such that the 
range of r is [to/2, +00). But thanks to the minimal character of the inner boundary 
S, we can extend (So, 7) to a larger Riemannian manifold (S ,7') with 7'| So = 7 
and 7' smooth at S. This is made possible by gluing a copy of S at S (cf. Fig. [3]). 
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Figure 4. Extended hypcrsurface £q depicted in the Kruskal-Szekeres 
representation of Schwarzschild spacetime. R stands for Schwarzschild radial 
coordinate and r for the isotropic radial coordinate. R = is the singularity and 
R = 2m the event horizon. Eq is nothing but a hypersurface i = const, where 
t is the Schwarzschild time coordinate. In this diagram, these hypersurfaces are 
straight lines and the Einstein-Rosen bridge S is reduced to a point. 



The topology of E is S 2 x R and the range of r in E is (0, +oo). The extended metric 
7' keeps exactly the same form as (|65|) : 

7^ dx l dx j = (l + y) 4 ( dr2 + r2de2 + r2 sin2 6d( f 2 ) ■ ( 67 ) 
By the change of variable 

2 

m 

r^r' = — (68) 
4r 

it is easily shown that the region r — > does not correspond to some "center" but is 
actually a second asymptotically flat region (the lower one in Fig. [3]). Moreover the 
transformation (|68|) . with 9 and ip kept fixed, is an isometry of 7'. It maps a point 
p of Eo to the point located at the vertical of p in Fig. [3] The minimal sphere S is 
invariant under this isometry. The region around S is called an Einstein- Rosen bridge. 
(Eg, 7') is still a slice of Schwarzschild spacetime. It connects two asymptotically flat 
regions without entering below the event horizon, as shown in the Kruskal-Szekeres 
diagram of Fig. 2J 



4-3. Bowen-York initial data 

Let us select the same simple free data as above, namely 

7« = fij, ixT = °> K = °> ^ = and = °- ( 69 ) 
For the hypersurface E , instead of R 3 minus a ball, we choose R 3 minus a point: 

E = M 3 \{0}. (70) 

The removed point O is called a puncture [3T] . The topology of E is § 2 x K; it differs 
from the topology considered in Sec. 14.11 (M 3 minus a ball); actually it is the same 
topology as that of the extended manifold E (cf. Fig. [3]). 
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Thanks to the choice the system to be solved is still (|4"4"l) - (l4"5l) . If we choose 
the trivial solution X = for Eq. (|4"5"1) , we are back to the slice of Schwarzschild 
spacetime considered in Sec. 14. 1[ except that now E is the extended manifold 
previously denoted Eg. 

Bowen and York [20] have obtained a simple non-trivial solution to the momentum 
constraint (|45j) (see also Ref. 15 ). Given a Cartesian coordinate system (x l ) = 
(x, y, z) on Eo (i.e. a coordinate system such that fij = diag(l, 1, 1)) with respect to 
which the coordinates of the puncture O are (0, 0, 0), this solution writes 

1 f y ijp i + ^) L ^\s ]X \ (7i) 



X' 



Ar 



where r := x 2 + y 2 + z 2 , e J k is the Levi-Civita alternating tensor associated with 
the flat metric / and (Pi,Sj) = (Pi, P2,Ps, Si, S2, S3) are six real numbers, which 
constitute the six parameters of the Bowen- York solution. Notice that since r / on 
£0, the Bowen- York solution is a regular and smooth solution on the entire So- 

The conformal traceless extrinsic curvature corresponding to the solution (f7Tj) is 
deduced from formula (|13p. which in the present case reduces to A lJ = (LX) 1 ^ ; one 
gets 



A l 3 — 

A "2r3 



'/" • ./••'/" - /'■' - 



P k x k 



1 (e* S*^' + ^Skx'j) , (72) 



where P l := f l ° Pj. The tensor A tJ given by Eq. (|72)l is called the Bowen- York extrinsic 
curvature. Notice that the Pi part of decays asymptotically as 0(r~ 2 ), whereas 
the Si part decays as 0(r~ 3 ). 

Remark : Actually the expression of A lJ given in the original Bowen-York article 
J20\/ contains an additional term with respect to Eq. \7$fy , but the role of this 
extra term is only to ensure that the solution is isometric through an inversion 
across some sphere. We are not interested by such a property here, so we have 
dropped this term. Therefore, strictly speaking, we should name expression {7%fy 
the simplified Bowen-York extrinsic curvature. 

The Bowen-York extrinsic curvature provides an analytical solution of the 
momentum constraint (|45[) but there remains to solve the Hamiltonian constraint 
(|44[) for ^ , with the asymptotic flatness boundary condition $ = 1 when r — > 00. 
Since I/O, Eq. (|4"4")l is no longer a simple Laplace equation, as in Sec. 14. 1[ but a 
non-linear elliptic equation. There is no hope to get any analytical solution and one 
must solve Eq. (j44[) numerically to get \& and reconstruct the full initial data (7, K) 
via Eqs. (l29l) - ([30"l) . 

The parameters Pi of the Bowen-York solution arc nothing but the three 
components of the ADM linear momentum of the hypersurface Eq Similarly, the 
parameters Si of the Bowen-York solution are nothing but the three components of 
the angular momentum of the hypersurface £0, the latter being defined relatively to 
the quasi-isotropic gauge, in the absence of any axial symmetry (see e.g. [51]). 

Remark : The Bowen-York solution with P l = and S l — reduces to the 
momentarily static solution found in Sec. \4-l\ *-e- is a slice t = const of the 
Schwarzschild spacetime (t being the Schwarzschild time coordinate). However 
Bowen- York initial data with P l = and S 1 7^ do not constitute a slice of Kerr 
spacetime. Indeed, it has been shown J^7f that there does not exist any foliation 
of Kerr spacetime by hypersurfaces which (i) are axisymmetric, (ii) smoothly 
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reduce in the non-rotating limit to the hyper surf aces of constant Schwarzschild 
time and (Hi) are conformally flat, i.e. have induced metric 7 = /, as the 
Bowen- York hyper surf aces have. This means that a Bowen- York solution with 
S' 1 =/= does represent initial data for a rotating black hole, but this black hole is 
not stationary: it is "surrounded" by gravitational radiation, as demonstrated by 
the time development of these initial data 122\ \4^j - 



5. Conformal thin sandwich method 



5.1. The original conformal thin sandwich method 

An alternative to the conformal transverse-traceless method for computing initial data 
has been introduced by York in 1999 |107j . The starting point is the identity 



(73) 



where N is the lapse function and (3 is the shift vector associated with some 3+1 
coordinates (f, x l ). The traceless part of Eq. (|73|) leads to 



A ij = — 
2N 

where A y is defined by Eq. 



>]). Noticing that 
2 ~ 



and introducing the short-hand notation 



7 



dt 



V 



we can rewrite Eq. (|?4"1) as 
1 



A 13 = 



2N 



1 



(74) 

(75) 
(76) 

(77) 
(78) 
(79) 
(80) 

Equation ([79"]) constitutes a decomposition of A^ alternative to the longitudi- 
nal/transverse decomposition (fT3"|) . Instead of expressing A 1 ^ in terms of a vector 
X and a TT tensor A^ T , it expresses it in terms of the shift vector (3, the time 

derivative of the conformal metric, 7 , and the conformal lapse N. 

The Hamiltonian constraint, written as the Lichnerowicz equation (|10|) . takes the 
same form as before: 



The relation between A ij and A 1 ? is [cf. Eqs. tfTJ-©] 

i« = >|/''.T'. 
Accordingly, Eq. ([77| yields 

2N L J 
where we have introduced the conformal lapse 



N := *" 6 iV. 



D,D*q 



R 



1 



-3 



^1 

12 



0. 



(81) 
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except that now A l i is to be understood as the combination (|79[) of j 3 and N. 
On the other side, the momentum constraint (fTTj) becomes, once expression (f?9"| is 
substituted for A l i , 

Dj (jWy 3 ) + Di - |* 6 ## - 16^. (82) 

In view of the system (l8TT) - (|82l) . the method to compute initial data consists in 
choosing freely 7y, 7* J , if, JV, £ and p l on S and solving (|5T|) - (f52"]) to get ^ and /3 Z . 
This method is called conformal thin sandwich (CTS), because one input is the time 

derivative 7 3 , which can be obtained from the value of the conformal metric on two 
neighbouring hypersurfaces S t and T, t +st ("thin sandwich" view point). 

Remark : The term "thin sandwich" originates from a previous method devised in 
the early sixties by Wheeler and his collaborators Q \101]j . Contrary to the 
methods exposed here, the thin sandwich method was not based on a conformal 
decomposition: it considered the constraint equations (QP-fiP as a system to be 
solved for the lapse N and the shift vector (3, given the metric 7 and its time 
derivative. The extrinsic curvature which appears in (Ip-flp was then considered 
as the function off, d~//dt, N and (3 given by Eq. \73\ ). However, this method 
does not work in general J3[. On the contrary the conformal thin sandwich method 
introduced by York \10T] and exposed above was shown to work [35]. 

As for the conformal transverse-traceless method treated in Sec. [3l on CMC 
hypersurfaces, Eq. (|82|) decouples from Eq. ([81~|) and becomes an elliptic linear equation 
for [3. 



5.2. Extended conformal thin sandwich method 



An input of the above method is the conformal lapse N . Considering the astrophysical 
problem stated in Sec. 12. 2[ it is not clear how to pick a relevant value for N. Instead 
of choosing an arbitrary value, Pfciffer and York [5D] have suggested to compute N 
from the Einstein equation giving the time derivative of the trace K of the extrinsic 
curvature, i.e. 

d_ 

dt 

3 



+ N 



D,D l N + 2D,ln^ D'N 



'( 

Att{E 



S) + A ij A i i 



(83) 



where S is the trace of the matter stress tensor as measured by the Eulerian observer: 
S = 7 A " y T A1I/ . This amounts to add this equation to the initial data system. More 



precisely, Pfciffer and York [80 suggested to combine Eq. (|83D with the Hamiltonian 
constraint to get an equation involving the quantity N 1 ^ — N 1 ^ 7 and containing no 
scalar products of gradients as the Di In ^D l N term in Eq. (|83|) . thanks to the identity 



D,D l N + 2D, In * D l N = 



(84) 



Expressing the left-hand side of the above equation in terms of Eq. (|83|) and 
substituting DiD 1 ^ in the right-hand side by its expression deduced from Eq. (fHTj) . 
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-R + ^T^ 2 * 4 + „i«i y *~ 8 + 2tt(E + 25)*" 
8 12 8 

>5 



we get 

Di&iNV 7 ) - (AT* 7 ) 

X - /3* DiJfJ * 5 = 0, (85) 
where we have used the short-hand notation 

and have set 

S := * 8 5. (87) 
Adding Eq. (f85|) to Eqs. (f8"Tj) and ([82]) . the initial data system becomes 



AD'f - -* + *~ 7 + 2t:E^- 3 - — * 5 = 

8 8 12 



Dj ( -(W ) + D j [jl - ^D L K = lfor? (89) 



f + -^ 2 * 4 + ^M '* -8 + 2tt(£ + 25)*" 
8 12 8 



A-D^iv* 7 ) - (iV* 7 ) 

K - f3 l A^srW 5 = 0, (90) 



where A^ is the function of N, [3\ % and ^ defined by Eq. (JTSJ). Equations p8 ]l -(|90 |l 
constitute the extended conformal thin sandwich (XCTS) system for the initial data 
problem. The free data are the conformal metric 7, its coordinate time derivative 7, 
the extrinsic curvature trace K, its coordinate time derivative K, and the rescaled 
matter variables E, 5 and p l . The constrained data are the conformal factor *, the 
conformal lapse N and the shift vector /3. 

Remark : The XCTS system f #<§ ) )- pffl ]) is a coupled system. Contrary to the CTT 
system t26]) - [WJ\ l, the assumption of constant mean curvature, and in particular 
of maximal slicing, does not allow to decouple it. 

5.3. XCTS at work: static black hole example 

Let us illustrate the extended conformal thin sandwich method on a simple example. 
Take for the hypersurface So the punctured manifold considered in Sec. 14.31 namely 

E =M 3 \{O}. (91) 

For the free data, let us perform the simplest choice: 

lij=fij, 7 J =0, K = 0, K = 0, E = 0, 5 = 0, and p l = 0, (92) 

i.e. we are searching for vacuum initial data on a maximal and conformally flat 
hypersurface with all the freely specifiable time derivatives set to zero. Thanks to 
(Pl). the XCTS system (JHHJ) - (EOJ) reduces to 

A* + -AijAv *~ 7 = (93) 
8 

%(i(£/3)^=0 (94) 

A(A>* 7 ) - -AijA^^N = 0. (95) 
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Aiming at finding the simplest solution, we notice that 

= (96) 
is a solution of Eq. ([M]). Together with f 3 = 0, it leads to [cf. Eq. {75]) ] 

A ij = 0. (97) 
The system (|93|) -([95 |) reduces then further: 

A* = (98) 

A(A* 7 ) = 0. (99) 

Hence we have only two Laplace equations to solve. Moreover Eq. ([98]) decouples 
from Eq. (|99|) . For simplicity, let us assume spherical symmetry around the puncture 
O. We introduce an adapted spherical coordinate system {x 1 ) = (r,9,tp) on So. The 
puncture O is then at r = 0. The simplest non-trivial solution of (|98|) which obeys 
the asymptotic flatness condition ^ — > 1 as r — > +00 is 

* = ! + (100) 
2r 

where as in Sec. I4.1[ the constant m is the ADM mass of So [cf. Eq. (|63|) ]. Notice 
that since r = is excluded from So, Vl/ is a perfectly regular solution on the entire 
manifold So. Let us recall that the Riemannian manifold (So, 7) corresponding to 
this value of \& via 7 = "F 4 / is the Riemannian manifold denoted (S ,"y) in Sec. 14.11 
and depicted in Fig. [3] In particular it has two asymptotically flat ends: r — > +00 
and r — » (the puncture). 

As for Eq. (f9"5| . the simplest solution of Eq. ([55]) obeying the asymptotic flatness 
requirement N^> 7 — > 1 as r — > +00 is 

jV* 7 = l + ~, (101) 
r 

where a is some constant. Let us determine a from the value of the lapse function at 
the second asymptotically flat end r — > 0. The lapse being related to N via Eq. (150)) . 
Eq. (jlOip is equivalent to 

/ fl\ , / a\ / tox -1 r + a 
N =i 1+ r)* =( 1+ r)( 1+ 2r) ^ 

Hence 

2a 

limAT=— . 103) 

r->0 m 

There are two natural choices for lim r ^o N. The first one is 

lim N = 1, (104) 

yielding a = m/2. Then, from Eq. (| 102[) AT = 1 everywhere on So- This value of N 
corresponds to a geodesic slicing. The second choice is 

lim N = — 1. (105) 

r— >0 

This choice is compatible with asymptotic flatness: it simply means that the 
coordinate time t is running "backward" near the asymptotic flat end r — > 0. This 
contradicts the assumption A > in the standard definition of the lapse function. 
However, we shall generalize here the definition of the lapse to allow for negative 
values: whereas the unit vector n is always future-oriented, the scalar field t is allowed 
to decrease towards the future. Such a situation has already been encountered for the 
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part of the slices t = const located on the left side of Fig. 2J Once reported into 
Eq. ([TD5]) . the choice (fTTiS]) yields a = -to/2, so that 



TON /. TO 



N =( i -fr)( i+ f r r 1 - (1 ° 6) 

Gathering relations (f9"u) . (|100p and (|106[) . we arrive at the following expression of the 
spacetime metric components: 

/l _ ™ \ 2 / mx 4 

g^dx^dx" = - (j-^r) dt 2 + [1 + —J [dr 2 + r 2 (dfl 2 + sin 2 0dy> )] . (107) 

We recognize the line element of Schwarzschild spacetime in isotropic coordinates. 
Hence we recover the same initial data as in Sec. 14.11 and depicted in Figs. [3] and [U 
The bonus is that we have the complete expression of the metric g on So, and not 
only the induced metric 7. 

Remark : The choices \104\ ) and U05\) for the asymptotic value of the lapse both lead 
to a momentarily static initial slice in Schwarzschild spacetime. The difference 
is that the time development corresponding to choice &104\ l (geodesic slicing) will 
depend on t, whereas the time development corresponding to choice H05\) will 
not, since in the latter case t coincides with the standard Schwarzschild time 
coordinate, which makes dt a Killing vector. 

5.4- Uniqueness of solutions 

Recently, Pfeiffcr and York [5T] have exhibited a choice of vacuum free data 

,K,K) for which the solution (^>,N,(3 l ) to the XCTS system (f8"8 > ([90" jl is not 
unique (actually two solutions are found). The conformal metric 7 is the flat metric 
plus a linearized quadrupolar gravitational wave, as obtained by Teukolsky [92], with 

a tunable amplitude. 7 3 corresponds to the time derivative of this wave, and both 
K and K are chosen to zero. On the contrary, for the same free data, with K = 
substituted by N = 1, Pfciffer and York have shown that the original conformal thin 
sandwich method as described in Sec. EH] leads to a unique solution (or no solution at 
all if the amplitude of the wave is two large). 

Baumgarte, O Murchadha and Pfciffer [TJ] have argued that the lack of uniqueness 
for the XCTS system may be due to the term 



7 , , _» 7 



-(N^)^A tJ A l ^- s ^-^ e ^ 7" + (W l kl + (W l (N* 7 )- 1 (108) 

o oZ L J L J 

in Eq. (|90p . Indeed, if we proceed as for the analysis of Lichncrowicz equation in 
Sec. 13.41 we notice that this term, with the minus sign and the negative power of 
(./V\E r7 )~ 1 , makes the linearization of Eq. (l90|) of the type DiD l e + ae = a, with a > 0. 
This "wrong" sign of a prevents the application of the maximum principle to guarantee 
the uniqueness of the solution. 

The non-uniqueness of solution of the XCTS system for certain choice of free data 
has been confirmed by Walsh [100] by means of bifurcation theory. 

5.5. Comparing CTT, CTS and XCTS 

The conformal transverse traceless (CTT) method exposed in Sec. |3] and the 
(extended) conformal thin sandwich (XCTS) method considered here differ by the 
choice of free data: whereas both methods use the conformal metric 7 and the trace 
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of the extrinsic curvature K as free data, CTT employs in addition A^ T , whereas 

for CTS (resp. XCTS) the additional free data is j J , as well as N (resp. K). 
Since Aj, T is directly related to the extrinsic curvature and the latter is linked to 
the canonical momentum of the gravitational field in the Hamiltonian formulation of 
general relativity, the CTT method can be considered as the approach to the initial 
data problem in the Hamiltonian representation. On the other side, 7 J being the 
"velocity" of 7 y , the (X)CTS method constitutes the approach in the Lagrangian 
representation [108J. 

Remark : The (X)CTS method assumes that the conformal metric is unimodular: 
det(7y) = / (since Eq. |7ff| j follows from this assumption), whereas the CTT 
method can be applied with any conformal metric. 

The advantage of CTT is that its mathematical theory is well developed, yielding 
existence and uniqueness theorems, at least for constant mean curvature (CMC) slices. 
The mathematical theory of CTS is very close to CTT. In particular, the momentum 
constraint decouples from the Hamiltonian constraint on CMC slices. On the contrary, 
XCTS has a much more involved mathematical structure. In particular the CMC 
condition does not yield to any decoupling. The advantage of XCTS is then to be 
better suited to the description of quasi-stationary spacetimes, since 7 3 = and 
K = are necessary conditions for dt to be a Killing vector. This makes XCTS 
the method to be used in order to prepare initial data in quasi-cquilibrium. For 
instance, it has been shown [57J S3] that XCTS yields orbiting binary black hole 
configurations in much better agreement with post-Newtonian computations than the 
CTT treatment based on a superposition of two Bowen-York solutions. Indeed, except 
when they are very close and about to merge, the orbits of binary black holes evolve 
very slowly, so that it is a very good approximation to consider that the system is in 
quasi-equilibrium. XCTS takes this fully into account, while CTT relies on a technical 
simplification (Bowen-York analytical solution of the momentum constraint) , with no 
direct relation to the quasi-equilibrium state. 

A detailed comparison of CTT and XCTS for a single spinning or boosted black 
hole has been performed by Laguna [68 . 

6. Initial data for binary systems 

A major topic of contemporary numerical relativity is the computation of the merger 
of a binary system of black holes [24] or neutron stars [84], for such systems are 
among the most promising sources of gravitational radiation for the interferometric 
detectors either groundbased (LIGO, VIRGO, GEO600, TAMA) or in space (LISA). 
The problem of preparing initial data for these systems has therefore received a lot of 
attention in the past decade. 

6.1. Helical symmetry 

Due to the gravitational-radiation reaction, a relativistic binary system has an inspiral 
motion, leading to the merger of the two components. However, when the two bodies 
are sufficiently far apart, one may approximate the spiraling orbits by closed ones. 
Moreover, it is well known that gravitational radiation circularizes the orbits very 
efficiently, at least for comparable mass systems [T5]. We may then consider that the 
motion is described by a sequence of closed circular orbits. 
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Figure 5. Action of the helical symmetry group, with Killing vector (.. Xt(P) 
is the displacement of the point P by the member of the symmetry group of 
parameter r. TV and (3 are respectively the lapse function and the shift vector 
associated with coordinates adapted to the symmetry, i.e. coordinates (t, x l ) such 
that d t = I. 



The geometrical translation of this physical assumption is that the spacetime 
(A4 , g) is endowed with some symmetry, called helical symmetry. Indeed exactly 
circular orbits imply the existence of a one-parameter symmetry group such that the 
associated Killing vector t obeys the following properties [46^ (i) I is timelike near 
the system, (ii) far from it, t is spacelike but there exists a smaller number T > such 
that the separation between any point P and its image Xt(-P) under the symmetry 
group is timclike (cf. Fig. [5]). £ is called a helical Killing vector, its field lines in a 
spacetime diagram being helices (cf. Fig. [5]). 

Helical symmetry is exact in theories of gravity where gravitational radiation does 
not exist, namely: 

• in Newtonian gravity, 

• in post-Newtonian gravity, up to the second order, 

• in the Isenberg- Wilson-Mathews (IWM) approximation to general relativity, 
based on the assumptions 7 = / and K = [6lJ 1102] . 

Moreover helical symmetry can be exact in full general relativity for a non- 
axisymmetric system (such as a binary) with standing gravitational waves [44] . But 
notice that a spacetime with helical symmetry and standing gravitational waves cannot 
be asymptotically flat [48] . 




To treat helically symmetric spacetimes, it is natural to choose coordinates (t, x l ) 
that are adapted to the symmetry, i.e. such that 




d t = e. 



(109) 



Then all the fields are independent of the coordinate t. In particular, 



7 J =0 and K = 0. 



(110) 
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If we employ the XCTS formalism to compute initial data, we therefore get some 
definite prescription for the free data 7 3 and K. On the contrary, the requirements 
piOp do not have any immediate translation in the CTT formalism. 

Remark : Helical symmetry can also be useful to treat binary black holes outside the 
scope of the 3+1 formalism, as shown by Klein \67^ , who developed a quotient 
space formalism to reduce the problem to a three dimensional SL(2, R)/SO(l, 1) 
sigma model. 

Taking into account () 1 1 0(1 and choosing maximal slicing (K = 0), the XCTS 
system ([SSJ-flSO) becomes 

D.D 1 ^ * + -AijA^ + 2ttE^- 3 = (111) 

8 8 

Dj (^(Lf3yA - 16nf = (112) 



N 

7\ / ATVTr7\ 



DiD l (N^') - (NV) 
where [cf. Eq. ([79])] 



J + lAijA***- 8 + 2ir(E + 25)*' 
8 8 



0, (113) 



A'* = -^r(L/3) ij . (114) 
2N 

6.2. Helical symmetry and IWM approximation 

If we choose, as part of the free data, the conformal metric to be flat, 

lii=fih ( 115 ) 
then the helically symmetric XCTS system (llll[) - (|113p reduces to 

A* + -A.jA 11 + 2itE^- 3 = (116) 



1, 
3' 

>7\ /J\7.Tr7\ 



A/3* + -V l Vjf3 j - {LffiVj In N = 16ttA^ (117) 



(118) 



where 



A ij = -L(Ll3) ij (119) 
2N 

and T> is the connection associated with the flat metric /, A := T>{D l is the flat 
Laplacian [Eq. gBJ)], and (L/3)« := V l ft +V j /3 i - §X> fe /3 fe /« [Eq. JISJ) with & = V% 
We remark that the system (jl 16[) - (II 18[) is identical to the system defining 
the Isenberg- Wilson-Mathews approximation to general relativity [6TJ 1102] (see e.g. 
Sec. 6.6 of Ref. [51]). This means that, within helical symmetry, the XCTS system 
with the choice K = and 7 = / is equivalent to the IWM system. 

Remark : Contrary to IWM, XCTS is not some approximation to general relativity: 
it provides exact initial data. The only thing that may be questioned is the 
astrophysical relevance of the XCTS data with 7 = /. 
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6. 3. Initial data for orbiting binary black holes 

The concept of helical symmetry for generating orbiting binary black hole initial data 
has been introduced in 2002 by Gourgoulhon, Grandclement and Bonazzola [521 157j . 
The system of equations that these authors have derived is equivalent to the XCTS 
system with 7 = /, their work being previous to the formulation of the XCTS method 
by Pfeiffer and York (2003) [80]. Since then other groups have combined XCTS with 
helical symmetry to compute binary black hole initial data [38l HI [2l E] . Since all 
these studies are using a flat conformal metric [choice (|115|) ]. the PDE system to be 
solved is (|116p - (lll8p . with the additional simplification E = and p 1 = (vacuum). 
The initial data manifold So is chosen to be R 3 minus two balls: 

So = M 3 \(£i U B 2 ). (120) 

In addition to the asymptotic flatness conditions, some boundary conditions must be 
provided on the surfaces Si and S2 of B\ and 82- One choose boundary conditions 
corresponding to a non- expanding horizon, since this concept characterizes black holes 
in equilibrium. We shall not detail these boundary conditions here; they can be found 
in Refs. [351 HOI SD (Ml ISS]. The condition of non-expanding horizon provides 3 
among the 5 required boundary conditions [for the 5 components ("J, N, The two 
remaining boundary conditions are given by (i) the choice of the foliation (choice of 
the value of N at S\ and S2) and (ii) the choice of the rotation state of each black 
hole ("individual spin"), as explained in Ref. [31 . 

Numerical codes for solving the above system have been constructed by 

• Grandclement, Gourgoulhon and Bonazzola (2002) [57] for corotating binary 
black holes; 

• Cook, Pfeiffer, Caudill and Grigsby (2004, 2006) [38j [31] for corotating and 
irrotational binary black holes; 

• Ansorg (2005, 2007) [TJH] for corotating binary black holes. 

Detailed comparisons with post-Newtonian initial data (either from the standard post- 
Newtonian formalism [T7] or from the Effective One-Body approach [2"5l 14"2"] ) have 
revealed a very good agreement, as shown in Refs. [4U1 131] . 

An alternative to (|120[) for the initial data manifold would be to consider the 
twice-punctured M 3 : 

So =K 3 \{0!,0 2 }, (121) 

where 0\ and O2 are two points of R 3 . This would constitute some extension to the 
two bodies case of the punctured initial data discussed in Sec. 15.31 However, as shown 
by Hannam, Evans, Cook and Baumgarte in 2003 [BO], it is not possible to find a 
solution of the helically symmetric XCTS system with a regular lapse in this cas^jjl 
For this reason, initial data based on the puncture manifold (|121|) are computed within 
the CTT framework discussed in Sec. [3] As already mentioned, there is no natural 
way to implement helical symmetry in this framework. One instead selects the free 
data Al^ T to vanish identically, as in the single black hole case treated in Sees. I4.ll and 
1431 Then 

= {LXfi. (122) 

I see however Ref. 1591 for some attempt to circumvent this 
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The vector X must obey Eq. (|45|) . which arises from the momentum constraint. Since 
this equation is linear, one may choose for X a linear superposition of two Bowen-York 
solutions fSec. 14. 3[) : 

X = X( P (i)_ S (i)) + X( P (2) S (2)j, (123) 

where X / P ( a ) jS ( a )\ (a = 1,2) is the Bowen-York solution (jTTj) centered on O a - This 
method has been first implemented by Baumgarte in 2000 [11] . It has been since then 
used by Baker, Campanelli, Lousto and Takashi (2002) [5] and Ansorg, Briigmann and 
Tichy (2004) [3J. The initial data hence obtained are closed from helically symmetric 
XCTS initial data at large separation but deviate significantly from them, as well 
as from post-Newtonian initial data, when the two black holes are very close. This 
means that the Bowen-York extrinsic curvature is bad for close binary systems in 
quasi-equilibrium (see discussion in Ref. |43j ) . 

Remark : Despite of this, CTT Bowen- York configurations have been used as initial 
data for the recent binary black hole inspiral and merger computations by Baker 
et al. JS| \99\j and Campanelli et al. \25\ \26\ \21\ 281. Fortunately, these initial 
data had a relative large separation, so that they differed only slightly from the 
helically symmetric XCTS ones. 

Instead of choosing somewhat arbitrarily the free data of the CTT and XCTS 
methods, notably setting 7 = /, one may deduce them from post-Newtonian results. 
This has been done for the binary black hole problem by Tichy, Briigmann, Campanelli 
and Diener (2003) [M], who have used the CTT method with the free data (7^ , A^ T ) 
given by the second order post-Newtonian (2PN) metric. This work has been improved 
recently by Kelly, Tichy, Campanelli and Whiting (2007) [66] . In the same spirit, 
Nissanke (2006) [75] has provided 2PN free data for both the CTT and XCTS methods. 

6.4- Initial data for orbiting binary neutron stars 

For computing initial data corresponding to orbiting binary neutron stars, one must 
solve equations for the fluid motion in addition to the Einstein constraints. Basically 
this amounts to solving V„T Ml/ = in the context of helical symmetry. One can then 
show that a first integral of motion exists in two cases: (i) the stars are corotating, 
i.e. the fluid 4- velocity is colincar to the helical Killing vector (rigid motion), (ii) the 
stars are irrotational, i.e. the fluid vorticity vanishes. The most straightforward way 
to get the first integral of motion is by means of the Carter-Lichnerowicz formulation 
of relativistic hydrodynamics, as shown in Sec. 7 of Ref. [50]. Other derivations have 
been obtained in 1998 by Teukolsky [S3J and Shibata |83] , 

From the astrophysical point of view, the irrotational motion is much more 
interesting than the corotating one, because the viscosity of neutron star matter is 
far too low to ensure the synchronization of the stellar spins with the orbital motion. 
On the other side, the irrotational state is a very good approximation for neutron 
stars that are not millisecond rotators. Indeed, for these stars the spin frequency is 
much lower than the orbital frequency at the late stages of the inspiral and thus can 
be neglected. 

The first initial data for binary neutron stars on circular orbits have been 
computed by Baumgarte, Cook, Scheel, Shapiro and Teukolsky in 1997 [El [13] in 
the corotating case, and by Bonazzola, Gourgoulhon and Marck in 1999 [19] in the 
irrotational case. These results were based on a polytropic equation of state. Since 
then configurations in the irrotational regime have been obtained 
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• for a polytropic equation of state [75J 023 IHZl [EH (HOI EE] (the configurations 
obtained in Ref. [5T] have been used as initial data by Shibata [53] to compute 
the merger of binary neutron stars); 

• for nuclear matter equations of state issued from recent nuclear physics 
computations [TBI [77] ; 

• for strange quark matter [751 172"] . 

All these computation are based on a flat conformal metric [choice (|115[) ]. by 
solving the helically symmetric XCTS system (|116[) - ([118[) . supplemented by an elliptic 
equation for the velocity potential. Only very recently, configurations based on a non 
flat conformal metric have been obtained by Uryu, Limousin, Friedman, Gourgoulhon 
and Shibata [98] . The conformal metric is then deduced from a waveless approximation 
developed by Shibata, Uryu and Friedman [55] and which goes beyond the IWM 
approximation. 

6. 5. Initial data for black hole - neutron star binaries 

Let us mention briefly that initial data for a mixed binary system, i.e. a system 
composed of a black hole and a neutron star, have been obtained very recently by 
Grandclement 55J and Taniguchi, Baumgarte, Faber and Shapiro [88, 89,. Codes 
aiming at computing such systems have also been presented by Ansorg [2] and Tsokaros 
and Uryu [95] . 
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